Delayed feedback as a means of control of noise-induced motion 
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Time-delayed feedback is exploited for controlling noise-induced motion in coherence resonance 
oscillators. Namely, under the proper choice of time delay, one can either increase or decrease 
the regularity of motion. It is shown that in an excitable system, delayed feedback can stabilize 
the frequency of oscillations against variation of noise strength. Also, for fixed noise intensity, the 
phenomenon of entrainment of the basic oscillation period by the delayed feedback occurs. This 
allows one to steer the timescales of noise-induced motion by changing the time delay. 

PACS numbers: 05.40.Ca Noise, 05.45.Gg Control of chaos, applications of chaos 
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Very often in practical application the need arises to 
control the properties of oscillations. Usually control as- 
sumes an enhancement in regularity of motion Q, Qi |^ . 
However, in some cases, for instance in medical applica- 
tions, one aims to disorder oscillations, since too strong 
coherence might be undesirable and even lead to dam- 
aging consequences, e.g. epilepsy or Parkinson's disease, 
0, 13 ■ During the last decade new methods for control of 
irregular self-sustained oscillations in deterministic sys- 
tems have been developed, including suppression of chaos 
by an external (periodic) signal [g, stabilization of un- 
stable periodic orbits embedded in a chaotic attractor 
by time-discrete control , or the use of a time-delayed 
feedback loop for the same purpose. 

Whereas the existing methods are designed to control 
deterministic oscillations or, most recently, noise-induced 
enhancement of deterministic oscillations [lol flll | and 
self-oscillations affected by noise ^2 ' there is a large class 
of systems that do not oscillate autonomously; but if they 
are forced even by a purely random process featuring no 
specific timescales, they demonstrate motion resembling 
a self-oscillatory process 0| . Prominent representa- 
tives of this class are excitable systems like neurons , 
chemical reaction systems |l6j| , and semiconductor nanos- 
tructures . The degree of closeness of their oscillations 
to ideally periodic ones, i.e. coherence, can depend reso- 
nantly on noise intensity Il3|. which is why it was called 
coherence resonance (CR) [ig- CR has been shown to 
occur in systems close to bifurcations J3 : in excitable 
systems [l^, and in bistable systems [l^. Remarkably, 
CR oscillators possess the fundamental property of self- 
oscillators, namely, the ability to synchronize [21.] . 

Frequently the timescale of oscillations in a CR system 
varies substantially depending on noise intensity. Since 
the latter is not easily controllable in practice, there is 
need to make a CR device robust against variation in the 
properties of noise. Another important task is to find 
a reliable way to deliberately change the timescales of 
noise-induced oscillations in a universal way without af- 
fecting intrinsic system parameters. Finally, the obvious 
need is to control the regularity of noise-induced motion. 
At present, all three problems remain a challenge. In the 
present Letter we propose to exploit time-delayed feed- 



back control to tackle all three issues. 

As a first example of a CR oscillator, we consider the 
noisy Van der Pol system closely before the Hopf bifur- 
cation, extended by a delayed feedback loop 
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(i^ - x^)y - Lulx + K{yr - y) + D^{t). 



Here, x and y denote phase variables at time t, while yr 
denotes the delayed variable y{t — T); K is the strength of 
delayed feedback, ^(i) is a random variable with Gaus- 
sian distribution, zero mean and unity variance, D is the 
noise intensity. We set the parameters v = —0.01 and 
ciJo = 1 at which a stable focus exists. 

First consider K = 0, i.e. no delayed feedback in Eq. 
(Q). While no hmit cycle occurs at £> = 0, apphcation 
of noise induces oscillatory motion as illustrated by the 
phase portrait in Fig. 1(a). The coherence of oscillations 
may be quantified by the correlation time tcor , estimated 
from the normalized autocorrelation function \l/(s) of y as 
tcor = /o°° |^'(s)|ds. In Fig. 2(a) the grey (green on-fine) 
line shows tcor vs D for K = Q. At small noise intensity 
D the oscillations are more coherent. 

In deterministic self-oscillatory systems, application of 
a delayed feedback in the form above acts as follows. If 
there exists an unstable periodic orbit of period T in the 
phase space, a feedback with delay time t — T can sta- 
bilize this orbit in some range of the control strength K. 
A CR system may have no periodic orbits, but noise may 
induce oscillations with a well-defined timescale that is 
associated with the spectral peak. We suppose that ap- 
plication of a delayed feedback can act by analogy with 
a system containing a periodic orbit, provided that r is 
equal or close e noug h to the basic period To of the noise- 
induced motion ^2] without feedback. Namely, it should 
suppress deviations from a reference state and thus en- 
hance the regularity of oscillations. 

To test this expectation, we switch on the control force 
in Eq. ((TJ. We set t = Tq with Tq = 6.17283951 « 
27r/wo- The phase portrait with control aX K = 0.2 is 
shown in Fig. 1(b) for the same D as in (a). Fig. 1(b) 
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FIG. 1: Phase portraits of noise-induced motion: (a),(b) Van 
der Pol oscillator at D = 0.003, (c),(d) FitzHugh-Nagumo sys- 
tem at D = 0.09 (the dashed lines denote the null- isoclines), 
(a),(c) if = 0; (b),(d) K = 0.2, r = To. 
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FIG. 2: (color on-line) Correlation time tcor vs. noise inten- 
sity D for (a) Van der Pol oscillator, (b) FitzHugh-Nagumo 
system. Grey (green on-line) lines: K = 0, black lines: 
K = 0.2, T = To- (b) grey (green on-line) circles: To for 
K ^0, black circles: Ti for K = 0.2, r = To. 



reveals a remarkable ordering of the oscillation as com- 
pared with Fig. 1(a). 

To quantify the ordering due to the feedback, we esti- 
mate tcor in dependence on D as above. It is given for 
K = 0.2 by the black line in Fig. 2(a). One can see that 
for any D the coherence of noise-induced oscillations be- 
comes larger when the delayed feedback loop is switched 
on and r is close to Tq. On the other hand, it was found 
that if T is far from an integer multiple of Tq, the coher- 
ence of oscillations, on the contrary, decreases. Also, for 
T = Tq the coherence of noise-induced oscillations in Eq. 

was found to increase monotonically with increasing 
K. In the following we fix if = 0.2. Next, we study how 
the feedback can affect the system's timescales. For this 
purpose, we set the noise at the intensity D — 0.003 as for 
Figs. 1(a), (b), and follow the evolution of Fourier power 
spectra with r, which is illustrated by Fig. 3(a). Without 
feedback (r = 0), system ^ has one pronounced peak /o 
in the spectrum. As r increases from zero, the peak fre- 



quency, height and width change. At t > 8, new peaks 
that change their positions, heights and widths with t 
become clearly visible. 

Since the control parameter of delayed feedback is the 
time interval t, we propose to describe the response of 
the system in terms of periods rather than frequencies. 
With feedback, we select all spectral peaks, and for each 
peak introduce the period T as the inverse of the peak 
frequency. Denote the period of highest peak by Ti. The 
dependences of T on r are given in Fig. 4 (a): Ti by cir- 
cles (yellow on-line), other T by crosses (blue on-line). 
One can see that variation of r changes Ti in a cer- 
tain range, doing so most effectively as r < To/2. As 
T increases beyond To/2, Ti drops quickly, and again 
increases with r with a similar slope as before. After 
T has increased by about Tq, Ti again drops abruptly 
to a lower branch, and again follows r, although with 
a smaller slope. These abrupt transitions to successive 
lower branches occur roughly every To time units, and 
each subsequent entrainment happens at a lower slope. 
The plot of Ti vs t exhibits a piecewise approximately 
linear dependence, the larger the r, the closer each seg- 
ment is to a straight line. 

The numerical results obtained above for the Van der 
Pol system can be understood in terms of a general the- 
ory of a canonical nonlinear oscillator with time-delayed 
feedback 

i + f{x,x) - K{ir - x) = 0. (2) 

Note that Eq. JQ) fits the form Eq. ^ if rewritten as 
a single second-order differential equation with f{x, x) = 



'{v — X )x + uj^x and y 



Without feedback (iC = 



or r = 0) the fixed point (xo, 0) is a stabe focus if 
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where partial derivatives are taken at the fixed point. In 
Eq. ^ xq = 0, and Eq. Q is true for the given param- 
eters. Setting K > does not change xq. Either with, 
or without feedback the noise-induced oscillations take 
place in the close vicinity of the fixed point. It is to be 
expected that the motion is influenced by the local prop- 
erties of this point. At r = the stable focus has a pair 
of complex conjugate eigenvalues Xq = po zL iq^, po < Oj 
qo / 0, and the value of go should give an estimate of the 
angular frequency. Indeed, the only peak of the power 
spectrum (Fig. 3(a)) has frequency /o ~ \qo\/2n. With 
T > 0, the system becomes infinite-dimensional, and pos- 
sesses a countable set of eigenvalues A. In order to ex- 
clude that the delayed feedback might induce the birth of 
a stable limit cycle via a Hopf bifurcation, thus provid- 
ing a trivial explanation for the remarkable ordering of 
oscillations, we perform a linear stability analysis of the 
fixed point of Eq. Following the standard routine 

of linearizing Eq. 121 around the fixed point the 
characteristic equation for A is derived: 
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Substituting X = p + iq, real and imaginary parts can be 
separated. The condition for a Hopf bifurcation is p — 0, 
q 0. Substituting it into the imaginary part of Eq. Q 
we obtain: 



cos qr 



K + df/dx 
K ■ 



(5) 



Since the right-hand side is larger than unity due to Eq. 
P|) , the Hopf bifurcation condition is not satisfied for any 
K and r. Thus, the delayed feedback in the given form 
cannot induce a Hopf bifurcation. 

The numerical solution of Eq. Q with f{x,x) from 
Eq. yields the eigenvalue spectrum X ~ p + iq a,s a 
function of r. The eigenperiods defined as T'^ — 27r/|g| 
(dots in Fig. 4(a)) coincide remarkably with the inverse 
peak frequencies of the power spectrum of the noise- 
induced oscillations as a function of t. The correspond- 
ing real parts p are given by dots in Fig. 4(b) (the seven 
largest p are shown). All p remain negative, but, as seen 
from Fig. 4(b), nonmonotonically change with r. As r 
increases, separate branches of p cross, thus providing a 
striking explanation of the strongly nonmonotonic, dis- 
continuous evolution of the dominant spectral peak of 
the noise-induced motion under delayed feedback: The 
period Ti of the highest spectral peak (circles (yellow 
on-line) in Fig. 4(a)) always coincides with the period 
T*^ of the least stable eigenmode, i.e. the one with the 
largest real part which we denote as pi (circles (yellow 
on-line) in Fig. 4(b)). The more stable eigenmodes result 
in the side peaks of the frequency spectrum. The more 
stable the modes are, the lower the peaks are. 

As pi oscillates with r, the degree of stability of the 
fixed point of the deterministic system is modulated, thus 
leading to modulation of the coherence of the stochastic 
motion, quantified by the correlation time tcor (solid line 
(green on-line) in the upper part of Fig. 4(b)). The local 
maxima of coherence occur when pi is close to zero, and 
Ti is close to Tq. 

The entrainment of Ti by t, which manifests itself in 
the almost piecewise linear dependence of Ti on r for 
large r, can be understood as follows. As shown above, 
it is related to the eigenvalue whose real part pi is closest 
to zero. Assuming pi w in Eq. we obtain Eq. jSJ. 
With df/dx = -t/ < if, and {K -u)/K k 1, this gives 
cos((7t) « 1 and \q\T « iim, where n is integer. Then 
the eigenperiod is 



(6) 



As illustrated by Fig. 4(b), pi is close enough to zero 
only for large r, for which the relation lO holds most 
accurately. To obtain the location of the maxima of pi , 
i.e. the maxima of coherence, substitute pi « 0, g = 
27m/r into the real part of Eq. Q), which yields r = 
2'Kn/LjQ — hTq. 

To summarize, delayed feedback applied to an oscil- 
latory system of the form gives rise to a countable 
set of eigenmodes of the fixed point, whose eigenperiods 




FIG. 3: Fourier power spectra of noise-induced oscillations in 
dependence on r for (a) Van der Pol oscillator at D = 0.003; 
(b) FitzHugh-Nagumo system at D = 0.09, K = 0.2. The 
spectrum is computed from y. 
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FIG. 4: (color on-line) Spectral peaks, coherence and eigen- 
values vs T a,t K = 0.2. (a),(b): Van der Pol oscillator at 
D = 0.003, (c),(d): FitzHugh-Nagumo system at D = 0.09. 
(a),(c) crosses (blue on-line): T, circles (yellow on-line): Ti, 
black dots: T*^. (b),(d) solid line (green on-line): tcor, (b) 
black dots: seven largest p, circles (yellow on-line): pi. 



and stability are controlled by r. The highest peak in 
the spectrum of the noise-induced motion is due to ex- 
citation of the least stable eigenmode. The coherence of 
oscillations is the higher, the less stable the mode is. The 
range of modulation of the peak frequency is largest for 
small tau and large K. 

Next, we consider another example of a CR oscillator, 
the FitzHugh-Nagumo system, which serves as a proto- 
type of an excitable system. Extending it again by a 
delayed feedback loop, we obtain 



dx 

'H 
dy 

dt 



= x + a + K{y,-y) + D^{t). 



(7) 



We set the parameters e — 0.01 and a = 1.1 such that 
a stable node is the only attractor of the system in the 
absence of feedback. 

Without feedback {K = 0), the mechanism for induc- 
ing oscillations by noise is different from the Van der 
Pol oscillator J^l- In Fig. 1(c), (d) dashed lines show 
the nuU-clines defined by dy/dt — (vertical) and by 
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dx/dt = (cubic parabola). They intersect at the fixed 
point, which is shghtly displaced to the left of the min- 
imum of the parabola for the parameters chosen. The 
nuU-clines divide the phase plane into four regions with 
different directions of phase velocity. Dots in Fig. 1(c) 
show the phase portrait with noise D = 0.09. 

In Fig. 2(b) the grey (green on-line) line shows tcor 
for Eq. |(7J) versus noise intensity, exhibiting a distinct 
maximum a.t D = 0.09. Also, grey (green on-line) circles 
show the basic period Tg of oscillations. Unlike Eq. 
here Tg changes substantially with noise, as was earlier 
shown in |20l |. 

Now, switch on the feedback with K — 0.2 and set t 
equal to the value of Tq = 4.12694 at optimum noise. 
The black line in Fig. 2(b) denotes tcor vs D, and shows 
that for any D the coherence of oscillations is higher with 
the feedback. However, this feature is not visible in the 
phase portrait (Fig. 1(d)). Black circles in Fig. 2(b) 
show the basic period Ti with feedback. It is evident 
that delayed feedback substantially reduces the variation 
of the noise-induced basic timescale. Note, however, that 
this may not be so if r is very different from Tg. 

Next, we study how the feedback can control the 
timescales and the regularity of noise-induced motion. 
Fix D at an optimum value 0.09, K at 0.2 and change 
T. The spectrum in dependence on r is given in Fig. 
3(b). With increasing r, the spectral peaks move to- 
wards zero, and the spectrum is gradually enriched by 
new peaks. As with the Van der Pol oscillator, select all 



visible peaks with periods T, and denote the period of the 
highest peak as Ti. In Fig. 4(c) T of several peaks are 
given by crosses (blue on-line), and Ti by circles (yellow 
on-line), depending on r. These dependencies are quali- 
tatively very much like those for Eq. |^ in Fig. 4(a). In 
Fig. 4(d) tcor is given depending on r, exhibiting oscilla- 
tory features. Local maxima of coherence occur when Ti 
is equal to Tq . Unlike in case of the Van der Pol oscilla- 
tor, the mechanism of delayed feedback control can not 
be explained by a local analysis of the fixed point since 
the oscillations are characterized by large excursions in 
phase space. Rather, a global analysis would be needed 
which is clearly beyond the scope of the present Letter. 

In conclusion, time-delayed feedback in the form of the 
difference between the current and a delayed state of the 
system can be used to control oscillations that are in- 
duced merely by noise. The most crucial parameter of 
such a control is the time delay r, depending on which 
the coherence of noise-induced oscillations increases or 
decreases. With this, a phenomenon of entrainment of 
the basic period of noise-induced motion by the time de- 
layed feedback is discovered. The latter ability is some- 
how reminiscent of classical synchronization phenomena, 
in that the externally imposed timescale r tunes the basic 
period of oscillations in the system, although it involves 
quite different mechanisms. 
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